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I. INTRODUCTION 


Extra dimensional brane-world models firstly introduced in [1] and then in [2] to resolve the 
hierarchy problem in fundamental interactions including gauge fields and gravity. The basic idea of 
these models and their followers is an assumption that matter and gauge fields are confined on a 3- 
dimensional brane embedded in a higher dimensional spacetime while the gravity, by definition, can 
travel in all dimensions. This is compatible with observations provided either the extra dimension 
volume to be in order of TeV scale as suggested in [1] or somehow warped over the brane as in [2]. 

The mentioned assumption is supported by string theory in which (D-)branes are defined to be 
where open strings ended. The latter correspond to the standard model fields in the low energy 
limit. In contrast, closed strings can propagate into extra dimension(s) off the brane and in the low 
energy limit correspond to gravitons. Despite of this justification, it is important to understand the 
mechanism of field localization on the brane in the low energy scales and independent of the string 
theory. This was not very satisfactory and all Standard Model’s particles could not be trapped 
into the brane in the original 5-dimensional version of this model, though there are successes by 
some 6-dimensional models [3] [4] [5]. 

In this regards, in [6], a 5-dimensional modihed Randall-Sundrum (MRS) model was proposed 
and shown that it improves the field localization behavior on the brane. The model looks like very 
similar to the original RS and even two geometries are locally the same, but it was shown that 
they differ globally [6]. An important difference is that, in contrast to the RS model, the new one 
does not contain any 5-dimensional cosmological constant. 

Regarding progress in the field localization by proposing [6], it is worth examining the effective 
4-dimensional gravitational behavior of this model. In this article, we investigate the gravitational 
perturbation of the background and find the graviton spectrum. We show that it excludes negative 
mass modes which signing stability of the model. Moreover its Newtonian gravity limit is studied. 

The organization of this paper is as follows. In section 2 we introduce the metric and remind 
its differences with the RS. In section 3, we perturb the background and find the graviton mass 
spectrum. In section 4, we study the Newtonian limit and conclude in section 5. 
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II. THE BACKGROUND METRIC 


Consider the 5-dimensional Einstein-Hilbert gravity action as in the following, 

^9 ~ 2 ^ j + Sbrane ( 2 - 1 ) 

in which S\,rane is a localized action of brane(s). Based on this action, the following metric can be 
introduced as a solution to the equations of motion which containes a 3-brane at the origin of fifth 
dimension [6]: 


+ dr^) (2.2) 

where r is the extra dimension varying from — oo to -|-oo and the Latin letters are from 1 to 5 
and Greek letters (brane coordinates) from 1 to 4. In addition; is Minkowskian metric with 
= diag{—l, 1,1,1). Henceforth, we use the nomination of [6] and call (2.2) the r-metric. 

It is worth recalling the well-known Randall-Sundrum (RS) metric as [7] 

dsjiS = + dz"^ (2.3) 

At the first sight it seems that the r-metric (2.2) can be converted to the RS metric (2.3), by the 
following coordinate transformation [6], 

g-fc|r| ^ l-k\z\ 
p-k\z\ 

dxf^ = - TT^dX^ (2.4) 

1 — k\z\ 

However, these two metric can not be the same spacetime for at least two reasons [6]: Firstly, 
the above coordinate transformation is singular at z = 1/k, so it can not be extended to all 
space. Secondly, a global coordinate transformation need to be an exact differential as dx^ = 
dx^/dX’^dX'^ + dx^/dzdz which is manifestly not the case. Moreover, it is shown in [6] that the 
consistent Einstein equation including the energy-momentum of the brane requires the cosmological 
constant to be zero as we already set in (2.1), in contrast to the Randall-Sundrum model which 
contains a negative 5-dimensional cosmological constant. Furthermore the field localization of 
different spin fields gets improvement compared to the original Randall-Sundrum model [6]. 
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III. FIELD EQUATIONS AND GRAVITON MODES 


To derive Einstein’s field equations for the matter of convenience, different coordinate system 
is adopted: 


ds^ = {rf^ydx^dx^) + dy^ (3.1) 

with a = ln(l — k\y\). By inserting (3.1) in the Einstein’s equation: 

= g^y{?>a + Qa‘^) 

= 6a^ 

where dot represents derivative with respect to the y extra dimension. 

To find graviton modes in the MRS model, metric would be purturbed around the fixed back¬ 
ground (3.1) as follows: 


(3.2) 

(3.3) 


Gab = GAB + hAB (3.4) 

Plug (3.4) into (2.1) and keep up to second order in h, then the variation due to , gives the 
linearized equation of motion [7]: 

NhMT - Nh + hrs - 

1 1 

+K‘^[TMTh% + TNTh%[ — -TjsANh — -QMNTTsh^^ - ^{‘^hMN “ fi'MTvh)] 

Kmn = 0 (3.5) 

where h = and Tmn is from (3.2) and (3.3). It could be shown that adopting the unitary gauge 
is feasible [8]: 

hf,5 = 0 k /i55 = /(x^)e”'^ := (p (3.6) 

where n is a constant. By rewriting the equation of motion (3.5) for different components we get: 
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/ii/ - component : -{dpd^h^u - d^d^hp^ - dyd^hpp + hpu) + Aa'^hp^ + -dpdji + -dpd^p 

+ — dpd^li — h — 4&h — dpO^p — Qa^h + 3&‘^p{2 + n)] 

3 3 ~ 

+ + 2hp^ - -Qpuh) - (3(7 + Qa‘^)hpy = 0 (3.7) 

1 ~ 3 

p5 — component ; -d^{dph — d^hp^) — -adpp = 0 (3-8) 

55 — component : -{d^d'^hpp — dpd^h) — -ah + (13 + 2n)&^p — 3a‘^h = 0, (3-9) 

in which h = g^'^hpy. Eq.s (3.7),(3.8) and (3.9) are a system of coupled partial differential equations 
which can be solved if one decouples them. To do this, we use the following tensor decomposition 
[9]: 


— ^Qpu + B (3.10) 

where E and B are scalars (called radion in the context of RS model), Cp is a divergenceless vector 
and Dpi, a divergenceless-traceless tensor (graviton). Then inserting in equations of motion one 
finds, 

pu — component : 2 ( ~ “^B^pSpv — Cp J’p — C^ p^p -\- Dp^^p -\- 2E pi, + {Egpi,)^^^ B^pi,^^ 

+ Cp ,u55 T Ci/^p^^ D-\~ ^p^pi/ -j- 4:a‘^ {—2Egpi, -\- B pi, + Cp^i, Ci,^p 

+ Dp,,) + ]^gp„ ( - 4E - - 4a{4E + Rf^^) - pP - 63^0“^ 

+ 3a‘^p{2 n)^ + a(^—gp„p — AgpyE + 2B py -|- 2Cp^i, + 2C„^p + 2,Dpp 

- \gpuB^p) = 0 (3.11) 

/i5 - component : - ^ae^^'f^p - ^Cp^u = 0 (3.12) 

55 — component : — -&g^'^B^p„ + (13 -|- 2n)a‘^e^^f — — 12a'^E — 6aE = 0 (3.13) 

One of the advantages of the decomposition (3.10) is making a distinction between the degrees of 
freedom of the model. Since we are interested in studying the effects of graviton, we can consistently 
set equal zero all parts of the decomposition (3.10) but the tensorial part Dpi,. Having done that. 
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the only non-trivial equation would be the — component equation of motion: 


dpd^D^y + Dpi, - ( 4 ct ^ + 2a)Dpy = 0 
Let us solve (3.14) by two boundary conditions. The first one is, 

lim Dp„{x,y) = 0 

and the second one is canceling out coefficients of delta Dirac in (3.14). 

By imposing the Fourier transform along the brane coordinates on eq. (3.14): 

+ — Dpi, + ( 4 ( 7 ^ + 2(7 )Dpy = 0 


(3.14) 


(3.15) 


(3.16) 


where p^ = g^'^PpPu = —m?. (3.16) is called Schrodinger-like equation and can be solved to find, 


Dpy - Wpyt 


+ U 




(3.17) 


where Wpy and Upy are some constant tensors and: 


i + yifFTi i-v/(fp+9 

"■ = - 2 - ’ = - 2 - 

with p^ := rjy^ppPy. Considering the second boundary condition results in: 


(3.18) 



ni) + (2 - n2 


= 0 


3 + a 

Wpy = - - Uay ■■= Fu 


3 — a 






(3.19) 


where 

The first boundary condition in which the solution (3.17) must vanish at y = ^ makes four 
divisions: 

1. = 0 then ni = 2, 77,2 = —1 and by condition (3.19), Upy = 0. So Wpy can be found by 


6 









normalizing the solution as follows, 


fi/fc 


\w 


flU \ 


i' 


f 

/ e = 1 


l-llk 

,(2ni-l)o- T _ 


dy = \ w 


flU 1 


‘ J (1 - k\y\y 


dy=l 


= 2 k 


For later use we find, 


liV 


= 2 k 


y=0 


(3.21) 


2. Either m? < 0 or 0 < < 2k‘^, then the solution (3.17) dose not meet the first boundary 

condition. So there is no graviton, he., = 0. 

3. 2k‘^ < m? < jk‘^, and by using (3.19), is proportional to so to find them, we 


normalize the solution to one; 




e ‘'Dl^dy = 1 


flU 


\ 


j (f2(1 - k\y\r + (1 - fe|2/|)-“ - 2F)dy = 1 


2(S + t^-2F 


Then 


fliy 


y=0 


= \u^,\\F-lf = k 


m? — 2k‘^ 
m? — k"^ 


(3.22) 


4. |A:^ < m? then n \^2 = (li*|o|)/2, using (3.19) and normalizing the solution to one, we find: 




e ""Dl^dy = 1 


/il/ 


\u 


yLU \ 


I (|F|2 + 1 - F(1 - A;|y|)+“ - F*(l - k\y\)-^‘^)dy = 1 


k 


2(|F|2 + l-2i?e( 


1+ia 
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Then 



+ 1 - 2fe(F)) = 


(3.23) 


Note that the final results in eq.s (3.22) and (3.23) are the same. 

Before going on it is worth to comment on these results. Firstly, there is no negative mass- 


squared graviton. This implies the stability of the background metric against tensorial perturba¬ 
tions. Our analysis reveals that there is a massless mode which as we will see in the next section is 


important to reach a Newtonian gravity limit. Moreover, there is a mass gap which separates the 
massless mode from the continuum. This is consistent with the fact that massive gravitons have 
not yet been found experimentally. 

IV. NEWTONIAN LIMIT 

The gravitational potential of a source mass m produced by the exchange of graviton reads as 
Yukawa potential [10]: 



(4.1) 


where sum is over the graviton spectrum found in the previous section and g is coupling constant 
for the gravitational interaction and can be given as [8], 



for m? = 0 
for > 2k‘^ 


(4.2) 


in which we have used (3.21)-(3.23). 

The interaction term in gravitational Lagrangian is [8]: 



(4.3) 


where is the energy-momentum tensor for matters on the 3-brane and such that: 




(4.4) 









because is contained in : 


e^'^^f\^\e--R^^Uydx + 


then: 


1 


1 


2K^k IGttG 


(4.5) 


where G is the 4-dimensional gravitational constant. Because of partly continuous mass spectrum, 
sum in (4.1) turns into an integration and it could be revealed that the integral measure dm is 
proportionate to k : 


Grn(2k ^ 

^ > 2k\r ^ r k 
runj^o 


Gm 


1 + 


/■+°° 1 m2 - 2k‘^ 
2A: m? -k^^ 


’’dm 


(4.6) 


the integration in right hand side is accounted as: 


r-+oo_2_2fc2 _ /2e’^^ - kEi{l, {V2 - l)kr)e^^W + kEi{l, {V2 + l)kr)e^^+^^^W 


'V2k 


m 

— k^ 


e-^'^dm = 


~ e 


-2y/2kr 


2{y/2-l)r V 

Then potential for r S> 1 at first orders reads as: 


Gm / ^- 2 ^/ 2 fcr 


-l-0(^)^ for fcr > 1 


(4.7) 


y(r-)~ _(l- 
r V 


r V 4(\/2-1)A;7' 


(4.8) 


So, as eq. (4.8) justifies, MRS model meets the gravitational Newtonian limit consistent with 
observations. 

V. DISCUSSION AND CONCLUSION 


In this paper, we constructed the gravitational perturbation of a MRS model (3.5) and found 
equations of motion for different components of metric perturbation as (3.7), (3.8) and (3.9). By 
the decomposition (3.10) and the unitary gauge, we solved the equations of motion. Mass spectrum 
for graviton modes depicted a massless mode which is in charge of the Newtonian gravity limit. 
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and a continuum that is responsible for small correction to the Newtonian limit in short distances. 
Interestingly, there is a mass gap which separates the single massless mode from the mass continuum 
that could be addressed why massive gravitons have not yet been detected experimentally. The 
higher dimensional corrections to the Newtonian gravity are exponentially suppressed due to the 
mass gap. Similar behavior has been recently reported in [11]. 

Finally it is worth mentioning that there is no negative mass-squared in the graviton spectrum 
which indicates the stability of the model. 
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